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Abstract 

The origin of magnetic flux noise in Superconducting Quantum Interference Devices with a power 
spectrum scaling as 1// (/is frequency) has been a puzzle for over 20 years. This noise limits the 
decoherence time of superconducting qubits. A consensus has emerged that the noise arises from 
fluctuating spins of localized electrons with an areal density of 5 x 10 17 m -2 . We show that, in 
the presence of potential disorder at the metal-insulator interface, some of the metal-induced gap 
states become localized and produce local moments. A modest level of disorder yields the observed 
areal density. 

PACS numbers: 03.67.Lx, 05.40. Ca, 73.20.Fz, 75.20.-g, 85.25.Dq 
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Well below 1 K, low-transition temperature Superconducting Quantum Interference De- 
vices [1] (SQUIDs) exhibit magnetic flux noise [2] with a temperature-independent spectral 
density scaling as l/f a , where / is frequency and 0.6 < a < 1. The noise magnitude, a 
few /i$oHz -1 / 2 at 1 Hz ($ is the flux quantum), scales slowly with the SQUID area, and 
does not depend significantly on the nature of the thin film superconductor or the substrate 
on which it is deposited. The substrate is typically silicon or sapphire, which are insulators 
at low temperature (T) [2J. Flux noise of similar magnitude is observed in flux and 
phase qubits. Flux noise limits the decoherence time of superconducting, flux sensitive 
qubits making scale-up for quantum computing problematic. The near-insensitivity of noise 
magnitude to device area j^, 0, f| suggests the origin of the noise is local. Koch et al. [3] 
proposed a model in which electrons hop stochastically between traps with different prefer- 
ential spin orientations. A broad distribution of time constants is necessary to produce a 1/f 
power spectrum 9]. They found that the major noise contribution arises from electrons 
above and below the superconducting loop of the SQUID or qubit [5|, 12], and that an areal 
density of about 5 x 10 17 m -2 unpaired spins is required to account for the observed noise 
magnitude. De Sousa [10j proposed that the noise arises from spin flips of paramagnetic 
dangling bonds at the Si-SiC>2 interface. Assuming an array of localized electrons, Faoro 
and Ioffe suggested that the noise results from electron spin diffusion. Sendelbach et 



al. [12J showed that thin-film SQUIDs are paramagnetic, with a Curie (1/T) susceptibility. 
Assuming the paramagnetic moments arise from localized electrons, they deduced an areal 
density of 5 x 10 17 m -2 . Subsequently, Bluhm et al. [3] used a scanning SQUID microscope 
to measure the low-T paramagnetic response of (nonsuperconducting) Au rings deposited 
on Si substrates, and reported an areal density of 4 x 10 17 m -2 for localized electrons. Para- 
magnetism was not observed on the bare Si substrate. 

In this Letter we propose that the local magnetic moments originate in metal-induced 
gap states (MIGS) [IjJ localized by potential disorder at the metal-insulator interface. At 
an ideal interface, MIGS are states in the band gap that are evanescent in the insulator and 

n 

extended in the metal 14j (Fig.l). In reality, at a nonepitaxial metal- insulator interface 
there are inevitably random fluctuations in the electronic potential. The MIGS are particu- 
larly sensitive to these potential fluctuations, and a significant fraction of them-with single 
occupancy-becomes strongly localized near the interface, producing the observed paramag- 
netic spins. Fluctuations 15( of these local moments yield T-independent 1/f flux noise. 
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To illustrate the effects of potential fluctuations on the MIGS we start with a tight-binding 
model for the metal-insulator interface, consisting of the (100) face of a simple-cubic metal 
epitaxially joined to the (100) face of an insulator in a CsCl structure (Fig. 2(a)). For the 
metal we assume a single s-orbital per unit cell and nearest neighbor (NN) hopping. For 
the insulator we place an s-orbital on each of the two basis sites of the CsCl structure and 
assume both NN and next-nearest neighbor (NNN) hopping. The parameters are chosen 
so that the metal s-orbitals are at zero energy and connected by a NN hopping energy of 
-0.83 eV. The onsite energy of the orbitals on the Cs and CI sites is taken to be -4 eV 
and 2 eV, respectively, and both the NN and NNN hopping energies are set to -0.5 eV. 
These parameters yield a band width of 10 eV for the metal, and 8 and 4 eV band widths, 
respectively, for the valence and conduction bands of the insulator with a band gap of 2 
eV (Fig. 2(d)). These band structure values are typical for conventional metals and for 
semiconductors and insulators. For the interface we take the hopping energy between the 
metallic and insulating atoms closest to the interface to be -0.67 eV, the arithmetic mean 
of -0.83 and -0.5 eV. 

The electronic structure of the ideal metal-insulator junction is calculated using a super- 
cell [l6] containing 20 x 20 x 20 metal unit cells and 20 x 20 x 20 insulator unit cells, a 
total of 24,000 atoms. The total density of states (DOS) of the supercell (Fig. 2(e)) shows a 
nearly flat DOS in the band gap region. The states in the insulator band gap are MIGS that 
are extended in the metal, decaying rapidly away from the interface into the insulator. Our 
model with a lattice constant of 0.15 nm yields an areal density of states for the MIGS of 
about 3 x 10 18 eV _1 m~ 2 , consistent with earlier self-consistent pseudopotential calculations 



To mimic the effects of interfacial randomness, we allow the onsite energy to fluctuate 
for both metal and insulator atoms near the interface [lsj]. Specifically we assume an energy 
distribution P(E) = (l/\/2ir5)exp[-{E - E ) 2 /25 2 }, where E is the original onsite energy 
without disorder, and 5 is the standard deviation. We characterize the degree of disorder by 
the dimensionless ratio R = 25 /W, where W is the bandwidth of the metal. For those MIGS 
that become localized, the energy cost, Ui, for double occupation is large, and we cannot 
use a noninteracting electron approach. Instead we adopt a strategy similar to that used by 
Anderson in his calculation of local moment formation 19( . We separate the space near the 



interface into 3 regions: (i) the perfect metal region (M), (ii) an interfacial region consisting 



3 



of 2 layers of metal unit cells and 2 layers of insulator unit cells (D) (Fig. 2(b)), and (iii) the 
perfect insulator region (I). Region (ii) is analogous to the impurity in Anderson's analysis. 

We first compute the single-particle eigenstates, of region D in isolation. For each 

of these states, we compute Uj (using a long-range Coulomb potential with an onsite cutoff 
of 10 eV) and the hybridization energy i -1 , due to hopping to the metal and the insulator {^(j] . 
With the computed values of Uj and r iy we solve Anderson's equation for the spin- dependent 
occupation for each localized state | i): 

i r Ep r 

(n ia ) = - dE'- ^- 5 . (1) 

Here, E ij(J = Ej i + Ui{rii a ) and a is the spin index. The net moment associated with the 
state is given by rrij = /xbK^o-) — (n%-a)\. Equation (DO) and the associated expression for 
the net moment of the localized states are calculated within the self-consistent Hartree-Fock 
approximation 19]. An m, ^ solution is obtained only when Uj/(Ep — Ej) exceeds a 
critical value which depends on FjjiEp — Ej). In the large Uj limit, it is more appropriate 
to start from the weak coupling limit {Tj = 0), where the localized state is populated by 
a single electron, and treat Tj as a perturbation. By calculating the areal density of such 
moment-bearing localized states we estimate the density of spin-| local moments. 

Figure 3 shows the calculated distribution p(E, U) in the isolated interfacial region for 
R— 0.05, 0.1, 0.15, 0.2, 0.25, and 0.3; for each value, higher values of U correspond to 
more localized states. As expected we see that, for any given degree of randomness, the 
states with energy inside the insulator band gap (the MIGS) or those at the band edges are 
most susceptible to localization. Figure 4 shows a perspective plot of the charge density of 
two states, with high and low values of Uj, showing the correlation between the degree of 
wavefunction localization and the value of Uj. Both states are centered in the insulator, a 
general characteristic of localized states in the band gap originating from the MIGS. 

Setting the Fermi energy at the insulator midgap value, we estimate the areal density of 
spins for a given degree of randomness R. The top panel in Fig. 5 depicts the distribution 
p(E, m) of the spin moments as a function of energy. We see that for small R virtually all the 
local moments are derived from the MIGS. The bottom panel of Fig. 5 shows the calculated 
areal density of local moments versus R. Our simple model thus indicates that moderate 
potential fluctuations (R ~ 0.15) at the interface produce an areal density of localized 



moments comparable to experimental values [2JJ. Although our analysis is for a specific 



model, we expect the general physical picture to remain valid for real materials. First, the 
formation of MIGS at a metal-insulator interface is universal, and their areal density is 

n 

rather insensitive to the nature of the materials as discussed in supplements 1201 and shown 

n 

numerically in Ref. [17| . Second, the formation of local moments from the combination of 
localized states and Coulomb interaction is a general phenomenon [ljj]. We also note that 
our analysis should not be significantly modified when the metal is superconducting. This 
is because the Ui for the localized states is generally much greater than the pairing gap. Of 
course, extended states with negligible U would be paired. 

Given our picture of the origin of the localized spin-| moments, how do they produce 
1/f flux noise with a spectral density S$(f) oc l// a ? The local moments interact via 



mechanisms such as direct superexchange and the RKKY interaction ll|, |22|, |23|, |24j between 



themselves, and Kondo exchange with the quasiparticles in the superconductor. This system 
can exhibit a spin-glass transition [25], which could account for the observed susceptibility 
cusp [12I near 55 mK. For T > 55 mK, however, experiments suggest that the spins are 



in thermal equilibrium |26j and exhibit a 1/T (Curie Law) static susceptibilit y [12L 1 13 ] - 
In this temperature regime, for hf « k B T standard linear response theory [23] shows 
that the imaginary part of the dynamical susceptibility T) = A(f,T)(hf/k B T). Here, 
A(f,T) oc Y l ^ a ,p p <x5(hf + E a - Ep)\((3 \ \ a)\ 2 , where is the /i-th component of 
the spin operator, a and (3 label the exact eigenstates, and P a is the Boltzmann distribution 
associated with state a. Combining the above result with the fluctuation-dissipation theorem 

n 

[151 ] which relates the flux noise to T), namely S<s>(f,T) oc (kBT/hf)x"(f,T), we 

conclude that the observed l/f a spectral density implies A(f,T) oc l// a (0.6 < a < 1). 
Assuming low frequency contributions dominate the Kramers-Kronig transform, this result 
is consistent with the observed 1/T static susceptibility, and the recent measurement 28 ] 
showing that flux noise in a SQUID is highly correlated with fluctuations in its inductance, 
However, without knowing the form of the interaction between the spins, one cannot derive 
this behavior for A(f, T) theoretically. 

In conclusion, we have presented a theory for the origin of the localized magnetic moments 
which have been shown experimentally to give rise to the ubiquitous low-T flux 1/f noise 
observed in SQUIDs and superconducting qubits. In particular we have shown that for 
a generic metal-insulator interface, disorder localizes a substantial fraction of the metal- 
induced gap states (MIGS), causing them to bear local moments. Although MIGS have 



5 



been known to exist at metal-insulator interfaces for three decades, we believe this is the 
first understanding of their nature in the presence of strong local correlation and disorder. 
Provided T is above any possible spin glass transition, experiments show that fluctuations 
of these local moments produce a paramagnetic x' an d a power-law, /-dependent x" which 
in turn leads to flux 1/f noise. It is important to realize that localized MIGS occur not 
only at the metal-substrate interface but also at the interface between the metal and the 
oxide that inevitably forms on the surface of superconducting films such as aluminum and 
niobium. There are a number of open problems, for example, the precise interaction between 
the local moments, its relation to the value of a, and the possibility of a spin glass phase 
at low temperature. A particularly intriguing experimental issue to address is why different 
metals and substrates evidently have such similar values of R, around 0.15. Experimentally, 
to improve the performance of SQUIDs and superconducting qubits we need to understand 
how to control and reduce the disorder at metal-insulator interfaces, for example, by growing 
the superconductor epitaxially on its substrate. 
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and S.G.L. thank M.Jain and J.D.Sau for fruitful discussions. This work was supported 
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Engineering Division, of the U.S. Department of Energy under Contract No. DE-AC02- 
05CH11231. S.C. acknowledges support from a Samsung Foundation. 
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Supplements 

Areal density of MIGS. We 



give a simple estimate of the areal density of MIGS. In 
a two-band tight-binding model [30], the amplitude squared of the evanescent solutions [31] 
close to the valence band edge has an energy- dependent decay length /3(E) = 2[2m*(E — 
Evbm)/^ 2 } 1 ^ 2 , where m* is the electron effective mass and E VBM is the energy of the valence 
band maximum. Near the conduction band edge Ecbm, P(E) = 2[2m*(EcBM ~ E)^ 2 ] 1 / 2 . 
The areal density N of MIGS in the insulator (in units of states per unit area) is given by 

Q 

/Ep poo fEp 1 

dE J dz V (E)e^ E > = 7] J dE— (2) 

where we have assumed the density of states r](E) of the metal to be constant over the 
energy range of the band gap. Inserting the expression for j3{E) into Eq.(j2J), we obtain 

N = V [(h 2 /2m*)(E F - E VBM )} 1 ' 2 (3) 



For most semiconductors and insulators |32j, m e /m* ~ 1 + C\/E g and E F — E VBM = C 2 E g 
with C\ ~ 10 eV and C 2 ~ 0.5; furthermore, for most metals r)(E) is of the same of order 
of magnitude. Consequently, the approximate expression 

NKTiiitfftmJdCz] 1 ' 2 (4) 

is relatively insensitive to the nature of both the metal and the insulator. Using the typical 
values 7}{E) w 2 x 10 28 m" 3 eV _1 and CiC 2 ~ 5eV, we obtain iV w 8 x 10 18 m~ 2 , in good 



agreement with pseudopotential calculation 17] for Al in contact with Si, GaAs or ZnS. 



Hubbard energy £/j. We calculate the Hubbard energy Ui for double occupation for 
states in the isolated D region by evaluating the integral 

Jd \r-r'\ 

over the supercell. Within our tight-binding supercell scheme, two additional factors need to 
be included, (i) The part of the Coulomb integral on the same atomic site is replaced with 
the value of an onsite Hubbard Uq. (ii) When the localization length (£) of the localized 
states is larger than the supercell size, there is overlap of wavefunctions from the neigh- 
boring supercell; this overestimates U{ for the very weakly localized states. Given that the 
participation number, Pj = 1/^,- ^(iv,-)! 4 ~ (^j/a) d in a disordered <i-dimensional system 
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with supercell lattice constant a and Ui oc 1/&, we map the U value of the finite supercell 



33| for f . 



onto that of an infinite supercell using a scaling law 

Hybridization energy broadening i^. The hybridization-energy broadening of the 
localized states arises from couplings to the extended states in the metal as well as those in 
the insulator, and is given by 

A = r? + r/ (6) 
r? = *\v % M LeP M {E), r} = niv/HJiE) (7) 

where p M ^\E) is the density of extended states in M (I) at the energy of the localized 
state E, and is the hopping matrix element between an extended state in M(I) and a 

localized state in D (ave indicates averaging over the extended states). Extended eigenstates 
in M(I) are a linear combination of constituent orbitals; the can then be expressed 

in terms of the coupling of these orbitals to those in D. For example, the localized states 
inside the band gap of the insulator are hybridized with only extended states in M, and 
Ti = Ef 1 w 7tV 2 di/W. (Here di is the charge of the localized state | i) in the unit cell layer 
immediately adjacent to M.) 
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FIG. 1: (Color online) (a) Schematic density of states, (b) MIGS at a perfect interface with energy 
in the band gap are extended in the metal and evanescent in the insulator. 
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FIG. 2: (Color online) (a) The metal (M) has a simple cubic structure with one atom per unit 
cell and the insulator (I) a CsCl structure with two atoms per unit cell, (b) Interfacial region (D) 
consists of 2 layers of metal unit cells and 2 layers of insulator unit cells. The lattice constant is 
0.15 nm. Computed DOS with Fermi energy (dotted red line) set to zero, (c) Typical metal with 
10 eV bandwidth, (d) Typical insulator with a 2 eV band gap separating two bands of about 8 
eV and 4 eV. (e) Metal-insulator interface with MIGS in the band gap of the insulator due to the 
presence of the metal. 
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FIG. 3: (Color online) Density of states distribution p(E, U) as a function of energy E and Hubbard 
energy U for 6 values of the randomness parameter R in the isolated D region of Fig. 2. For a 
given value of R, the highest values of U, resulting in the most highly localized states, appear in 
the band gap of the insulator and at the band edges. The position of the insulator band gap is 
represented by black dashed lines. 
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FIG. 4: (Color online) Perspective view images of the 2-dimensional probability density distribu- 
tion at the interfacial region (D) along directions parallel to the interface (x- and y-directions), 
integrated along the z-direction. (a) States with [7j=3.25 eV and 2£j=-0.24eV and (b) with C/j=0.35 
eV and i?j=-0.23eV, respectively. 
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FIG. 5: (Color online) (a) Electron density distribution p(E,m) for 6 values of R. We simulated 
5000 different configurations of disorder for each value of R. The position of the insulator band 
gap is represented by black dashed lines. Virtually all the magnetic moments are from the MIGS 
in the band gap of the insulator, (b) Integrated spin density versus randomness parameter R. For 
R = 0.05, we estimate the spin density to be less than 0.01 x 10 17 m~ 2 . 



13 



